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Abstract 

The spontaneous breakdown of SO(IO) symmetry of the IIB matrix model has been studied by using the 
improved mean field approximation (IMFA). In this report, the eighth-order contribution to the improved 
perturbative series is obtained, which involves evaluation of 20410 planar two-particle irreducible vacuum 
diagrams. We consider SO(J)-preserving configurations as ansatz (d = 4, 7). The development of plateau, 
the solution of self-consistency condition, is seen in both ansatz. The large ratio of the space-time extent of 
fif-dimensional part against the remaining (10 - (i)-dimensional part is obtained for S0(4) ansatz evaluated 
at the representative points of the plateau. It would be interpreted as the emergence of four-dimensional 
space-time in the IIB matrix model. 
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I. INTRODUCTION 



Superstring theory is supposed to provide a unified microscopic description of the universe 
including gravity. Recent progress revealed that many perturbative vacua thus far examined are 
related to each other, which enforces us to pursue non-perturbative formulation of superstring 
theory. One of those constructive definitions is the IIB matrix model |ll] formulated in a form of 
the large-A'^ reduced model of the ten-dimensional supersymmetric SU(A^) Yang-Mills theory. A 
characteristic feature of this model is that the eigenvalue distribution of ten bosonic matrices is 
interpreted as space-time, i.e., the space-time itself is treated as a dynamical variable of the model. 
It opens us a way toward explaining the origin of our four-dimensional space-time in the context 
of superstrings ibt. 

In this report, we consider the possibility of spontaneous breakdown of the original SO(IO) 
symmetry through the dynamics of the IIB matrix model. It is examined by the anisotropy of 
eigenvalue distributions of the matrices. In this regard, the extents of space-time measured by the 
second moment of the matrices vary according to directions. It may occur that the J-dimensional 
subgroup of the rotational symmetry stays intact, which reflects that J-dimensional space-time 
emerges as the vacuum of the IIB matrix model. 

Exploring non-perturbative dynamics of a model is not an easy task in general. Monte Carlo 
simulation, a powerful tool that has been successfully applied to various physical models, is 
plagued by the sign problem in this case, because the action has a complex phase originated from 
the fermionic part.' Whilst it is suggested that the complex phase plays crucial role in the sponta- 
neous breakdown of Lorentz symmetry Isl. 

Here, we instead exploit a technique called the improved mean field approximation (IMFA). It 
is one of variational methods capable of exploring non-perturbative solutions of the model. It was 

n n 

developed in Ref. j5<\, and first applied to IIB matrix model in Ref. la], in which a systematic im- 
provement to higher order terms was proposed and evaluated up to third order. The mechanism of 
spontaneous breakdown of rotational symmetry was further examined with the simplified models 



The analysis of the IIB matrix model by the IMFA method was proceeded to incorporate even 
higher-order contributions, up to fifth order in Ref. ||8|] and up to seventh order in Ref. ||9|]. In the 
latter work the automated calculation procedure was also developed. 



^ A new technique called the factorization method is proposed to resolve the complex action problem in the Monte 
Carlo simulations 
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The series of works showed that the spontaneous breakdown of Lorentz symmetry would occur 
to resuk in the emergence of the four-dimensional space-time. It was also observed that the extent 
of space-time of the four-dimensional part is larger than that of the remaining six-dimensional part, 
which is suggestive of the compactification scenario of string theory realized within the context of 
non-perturbative dynamics of the IIB matrix model. 

In the IMFA prescription, the Gaussian terms with auxiliary parameters are introduced into the 
action of the original model, which may be considered as mean fields in the mean field approxi- 
mation. The self-consistency condition for those parameters is given by the principle of minimal 
sensitivity |lfl|l as a guide line. It should be realized as a region of parameter space denoted as 
plateau. The development of plateau provides the indication whether or not the approximation is 
working well. 

In the present report, we proceed with the IMFA analysis of the IIB matrix model up to the 
eighth-order contribution. Since we are particularly interested in the spontaneous breakdown of 
SO(IO) symmetry, we restrict ourselves to two particular choices of breaking patterns, namely, 
S0(4) ansatz and S0(7) ansatz, in which four- (seven-) dimensional rotational symmetry stays 
intact. By incorporating higher-order terms we expect a clearer signal for the emergence of plateau, 
as well as better estimates of the free energy and the ratio of the space-time extent in the true 
vacuum of the model. 

The paper is organized as follows. In Section|lll we provide a description of the model and the 
application of the IMFA analysis in a concrete procedure. The particular choices of configurations 
to be examined are given in Section|lIll In Section|lVlwe present the results. Section|V|is devoted 
to the conclusion and discussions. 

II. IIB MATRIX MODEL AND IMFA ANALYSIS 
A. IIB Matrix Model 

The model we are considering here is the IIB matrix model defined by the partition function. 



where bosonic variables and fermionic variables if/a are both traceless N x N Hermite matri- 
ces. Qi = 1, . . . , 10) transforms as a vector under SO(IO) rotation, while i/^q, (a = 1, . . . , 16) 




(1) 
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transforms as a left-handed spinor. The action has symmetry under the U(A^) matrix rotation, 
ten-dimensional Lorentz symmetry, and the type IIB supersymmetry. 

The coupling constant go may be absorbed by the rescaling of the fields. Then the action takes 
the following form, 

S = A^Tr|-i[A^,A,]' - i^r'^[A^,^]| , (2) 

where we have chosen ^q^A^ = 1, and and denote the rescaled fields. We will examine the 
large- behaviour of the model. 

The extent of space-time along a certain direction ju is defined by the second moment of the 
matrix: 

^/ = ^{TrA,^). (3) 
B. Improved Mean Field Approximation 

The action of the IIB matrix model Q does not have quadratic terms, and therefore the ordinary 
perturbation theory is not directly applicable. To such cases we employ a technique called the 
improved mean field approximation (IMFA) that is elucidated in the following. 

We introduce a quadratic term 5o with arbitrary parameters collectively denoted by mo and 
deform the original action as 

S So + A(S-So), (4) 

where Aha nominal parameter that should be taken to 1 . The second term of the deformed action 
© may be considered as an interaction term with a coupling constant A. So, we can formulate a 
perturbation theory and obtain a power series expansion with reference to A. The series up to some 
finite order n would provide an nth order approximation of the original model at A = 1 . By this 
procedure, a formal perturbative expansion of a model can be constructed even when the model 
does not have quadratic terms, for example, as in the case of the IIB matrix model. 

The perturbative series thus obtained depends on the arbitrary parameters niQ. We have to 
determine the values of those parameters by some means. We adopt here the principle of minimal 
sensitivity 1 1^ as a guiding principle; the true value of a physical quantity should be given when it 
depends least on the arbitrary parameters. It is because the deformation © of the action becomes 
trivial by construction if A is taken to 1 , and the result is independent of those artificially introduced 
parameters. The dependence is brought in due to the truncation at the finite order of perturbation. 
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Therefore, if a region of parameter space exists in which the series becomes almost constant, 
there the parameter dependence vanishes effectively, and the true value should be reproduced. 
We call this region "plateau". It has been tested and shown to work on a number of systems 




\11 



m. 



The concept of plateau above is rather obscure, and thus we need a more concrete criterion for 
the distinction of plateau in such a manner that the ambiguity due to subjectivity of recognition 
be excluded as much as possible. It is usually seen that the series as a function of the arbitrary 
parameters fluctuates on and near the plateau about the true value, and it is accompanied by a 
number of extrema of the function. This leads to a practical criterion for identifying plateau by the 
accumulation of extrema of the series with reference to the parameters.^ It has been adopted in the 
previous works The values of the series are estimated at the extrema as the representatives 

of the plateau, which should give good approximations if the series is convergent. We adopt this 
criterion in the present analysis. 



C. Application to IIB Matrix Model 

Let us proceed to the IIB matrix model and evaluate the free energy of the model by the IMFA 
method. We introduce a quadratic term 5 o of most generic form that preserves U(A^) symmetry, 

5o = NTrl^^M^A,A, + ^m^yp^Y'^^'ip^ . (5) 

Here, M^y and m^^p are arbitrary parameters. M^y are symmetric with /i and v, while m^yp are 
totally anti-symmetric with /i, v and p. 

The prescription of the IMFA method is alternatively formulated by starting with the "massive" 
theory defined by the action, 

5' = 5o(M,m) + S(i), (6) 
where S {X) is given by inserting a formal parameter X to the original action Q as 

S = A^Tr|-^[A„Ay]' - J^^rnA„<A]| . (7) 

^ An alternative approach for the identification of plateau from the profile of functions with the help of histograms is 
proposed in Ref. L121 . 
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A physical quantity / is evaluated by the perturbative series expansion in terms of the coupling 
constant A up to some finite order n: 



mM,m) = Y,A'UM,m). 



(8) 



*:=0 



Then, we perform the following substitution of parameters, 



M. 



(l-^)M^v, 

m^^vp — > (1 - A)mf,yp , (9) 

and rearrange the series in powers of A. Finally, We disregard OiA"""^) terms and set the formal 
expansion parameter Aio \. Following those steps, we obtain the improved series f{A; m). It will 
be denoted as 



/ f(M, m) = / U; (1 - A)M, (1 - A)m) . (10) 

As is seen in Eq. ®, the model is shifted to "massless" case at A = 1, and thus the original model 
would be reproduced. 

For a technical reason, we exploit the fact that the ordinary free energy is related to the two- 



particle irreducible (2PI) free energy by the 



egendre transformation, the latter of which is ex- 



pressed by the sum of 2PI vacuum diagrams lll3ll . It is because the number of diagrams incorpo- 
rated reduces drastically. 

We evaluate the sum of 2PI vacuum diagrams in terms of the exact propagators. Since we are 
interested in the large- A'^ limit, the diagrams to be evaluated are restricted to planar ones. The 2PI 
free energy G is given by: 

t2 



G(C, u)/N' = 3(1 + log 2) + <{ -^tr log C + ^trlog 



+^|-i(tr(c2) - (trC)') - ^c,Mir^inj + ..., 

where C^y and = Uf^ypY^'^P /3\ are the exact propagators of and i/r, respectively: 

((^A')'/^v)/t/) = —C^y6ii6jk , 

1 i 



(11) 

(12) 
(13) 



C is the charge conjugation matrix. The additive constants are adjusted according to the definition 
in Ref. We then perform the Legendre transformation to obtain the free energy F by 



F(M, m) 



G(C, + Mf,yCf,y m^vpUf^yp 



(14) 



C=C(M,m),u=u{M,m) 



where C{M, m) and u{M, m) are determined by the solutions of the following relations: 

dG{C,u) 1 ^ ^ dG{C,u) 1 ^ 

+ xM^v = 0, — -m^vp = 0. (15) 

The improved free energy F is obtained by applying the procedure dTUI) to Then, we search 
for the extrema of F{M, m) with reference to the parameters M^y and m^^p, and identify the plateau 
by the accumulation of them. The values of free energy and other physical quantities are evaluated 
at the extrema. 



III. ANSATZ 



In the case of the IIB matrix model, a huge number of parameters are introduced along with 
the quadratic terms: 10 real numbers for M^y (assumed to be diagonalized), and 120 for m^yp. It 
may demand enormous efforts to explore the whole parameter space for the solutions of plateau 
condition. Instead, we restrict ourselves to particular configurations in which SO(J) subgroup of 
SO(IO) stays intact. The explicit forms of configurations are chosen according to the guideline 
described in Ref. [8]. In this report we concentrate on J = 4 and d = 1 cases in particular. 

Since we evaluate 2PI free energy in the first step, we specify the forms of exact propagators in 
each ansatz- 

SO(7) ansatz. The fermionic propagators are represented by the rank three anti-symmetric tensor 
Mpyp. A non-zero element of u accompanies three-dimensional block by considering the permuta- 
tion of indices, and thus the SO(IO) symmetry breaks down to S0(7) x S0(3). 

In this case, the propagators take the following forms with three parameters, ci, C2, and u: 



ili = uT'''''^'\ (16) 



SO(4) ansatz. We assume that S0(4) symmetry is preserved. The remaining six-dimensional 
subspace is decomposed into two three-dimensional blocks. This leads to the pattern of breaking 
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FIG. 1: Distributions of extrema of the improved free energy in the parameter space for S0(4) ansatz (left) 



and S0(7) ansatz (right). 



1 /m 




1/m 



1/mi 




'■2 1 / 

1 l/"i2 



1/mi 



as S0(4) X S0(3) x S0(3) x Z^. The extra Z2 factor derives from the permutation among two S0(3) 
factors with the reversion of the first direction (to preserve parity in total). 

In this case, the propagators are represented by the form given below with three parameters, ci, 
C2, and u: 




IV. RESULTS 

We evaluated the free energy of the IIB matrix model for S0(4) ansatz and S0(7) ansatz by 
the IMFA method up to eighth order of series expansion. The number of planar 2PI diagrams are 
shown in TableU In this report, the eighth-order contribution consisting of 20410 distinct diagrams 
are newly evaluated. The automated procedure to generate and evaluate the set of diagrams has 
been developed in Ref. lO], which is used in the present study with slight extensions. 

The improved free energy is obtained by following the procedure described in Sec. HH through 
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FIG. 2: Free energy evaluated at the extrema for S0(4) ansatz (left) and S0(7) ansatz (right). Horizontal 
axis represents the order of approximation. 



8 
S 
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S0(7) ansatz 



FIG. 3: Ratio of the extent of space-time evaluated at the extrema of the improved free energy for S0(4) 
ansatz (left) and S0(7) ansatz (right). Horizontal axis represents the order of approximation. 
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S0(7) ansatz 



the 2PI free energy G. Then, we search for the extrema of the improved free energy in the space 
of parameters mi, mj, and m (dual to ci, C2, and u, respectively). The extrema of free energy for 
the eighth-order improved series are listed in Table Ull We also confirmed that all the extrema of 

nil 

lower orders found in the former works BalSD are reproduced except one in seventh order of S0(4) 
ansatz. 

The distributions of extrema in the parameter space are plotted in Fig. [T] It is observed that a 
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TABLE I: Number of planar 2PI vacuum diagrams of ("massive") JIB matrix model. 



Order 


Number of diagrams 


0th 


2 


1st 


2 


2nd 


2 


3rd 


4 


4th 


12 


5th 


49 


6th 


321 


7th 


2346 


8th 


20410 



small region exists for both ansatz in which a number of extrema are found close to each other 
(shown in a dashed circle). It is recognized that extrema of different orders as well as those of 
respective orders belong to this accumulation. There are also extrema forming a line which would 
be considered as overshoots that characteristically appear around the edge of a plateau. 

In Fig. 121 the values of the free energy evaluated at the extrema are plotted for S0(4) and S0(7) 
ansatz against the order of the IMFA analysis. The bullets (•) correspond to the extrema that 
belong to the accumulation, while the circles (o) correspond to the other extrema. The cross marks 
(x) represent the unphysical extrema at which ci < or C2 < 0. 

Next, we consider what shapes of configurations are realized as the solutions. It is examined 
by evaluating the extent of space-time defined as the moments of the eigenvalue distributions: 

dF 



= i{TrA,A,) 



I dm 



dF 



improved 

(18) 



improved 



Here, the notation . ^ denotes the application of the IMFA procedure to the series on the left. 

Improved ^ ^ ^ 

The values of R, r and the ratio p = R/r evaluated at the extrema of free energy are also listed 
in Table im for each ansatz. The ratio is plotted in Fig.|3]for the order of approximation (the marks 
are the same as those of free energy). They show distinctive behaviour between S0(4) and S0(7) 
ansatz. In S0(7) case the ratio stays at p ^ 2.5, while in S0(4) case it grows larger as the order 
increases. It should be remarked that we obtained the ratio of the space-time extent in S0(4) ansatz 
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TABLE 11: Numerical values. 



Ansatz 


Free energy 






P 


S0(4) 


2.646862 


11.400026 


0.033561 


18.430393 




-4.056151 


11.234590 


0.035111 


17.887798 




-2.393029 


6.430896 


0.078981 


9.023488 




6.507018 


6.249471 


0.171968 


6.028345 




-29774.524354 


165.895937 


8.280539 


4.475985 




-3610.819690 


23.506193 


532.358893 


0.210130 


SO(7) 


2.885793 


0.831046 


0.123696 


2.591994 




2.931477 


0.898519 


0.148787 


2.457429 




-52.619043 


1.116045 


9.571141 


0.341475 




-108.828989 


0.397967 


6.643743 


0.244747 



at an order of magnitude larger than the isotropic configuration, which has not been seen until the 
present eighth-order calculation. 



V. CONCLUSION AND DISCUSSIONS 



We applied the improved mean field approximation (IMFA) to the HB matrix model and eval- 
uated up to eighth order contribution for configurations preserving SO(J) rotational symmetry as 
ansatz. We examined d = A and d = 1 cases in the present study. 

In order to solve the consistency condition for the arbitrary parameters called plateau condition, 
the stationary points (extrema) of the improved free energy are searched. It is observed that a 
region exists for each ansatz in which a number of extrema of eighth order as well as the extrema 
of different orders gather close to each other. It may be considered as a development of plateau. 
In the former works the discrepant behaviour was seen between even orders and odd orders of 
the improved series. There were no extrema found for S0(7) ansatz at even orders, and the free 
energy at extrema for SO(4) ansatz gave somewhat different sequence of values. Such discrepancy 
tends to be resolved with eighth-order terms taken into account. It would be consistent with the 
speculation that the result should be irrespective of the order of approximation, if the series is 
convergent and high enough orders are incorporated. 
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We evaluated the extent of space-time of J-dimensional part and that of (10 - J)-dimensional 
part in SO(d) ansatz. They are given by the moments of the eigenvalue distributions at the extrema 
of free energy as the representative points of plateau. The deviation of the ratio of the extents from 
1 represents the anisotropy of space-time realized as a non-perturbative vacuum of the IIB matrix 
model, and thus provides an indication of the spontaneous breakdown of Lorentz symmetry. In 
the present work, the ratio of an order of magnitude larger than the isotropic configuration was 
obtained for SO(4) ansatz. It might suggest that the actual ratio realized in the true vacuum of 
the theory may possibly be infinite, which implies the emergence of extended four-dimensional 
universe with the remaining six-dimensional part being compactified, as seen in the universe. 

Still we do not have clear signal enough to infer the formation of plateau, and therefore no 
reasonable estimate for free energy nor other physical values would be assured at this stage. The 
configuration with smaller values of free energy should give dominant contribution in the whole 
configuration space, and it is supposed to be realized as our universe. In this regard, we can not yet 
definitely tell which of the ansatz would be plausible. More sophisticated scheme for identifying 
plateau and extracting the physical quantities would be required along with the effort to proceed 
to even higher order contributions. 
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